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D ■ 

|Jh I Abstract. We discuss numerical results for the 1 — D spin 1/2 antiferromagnetic 

#>f-\ ■ Heisenberg model with next-to-nearest neighbour coupling and in the presence of an 

uniform magnetic field. The model develops zero frequency excitations at field depen- 
^-H , dent soft mode momenta. We compute critical quantities from finite size dependence 

^ ■ of static structure factors. 

S ' 1 Introduction 

CN ■ 

O; 

t^^ I The study of quantum spin systems in one space dimension (spin chain) can 

0^ ■ be motivated as follows: There are models which are exactly soluble, e.g., the 

-M I Heisenberg model. "Exact" numerical results can be obtained by diagonalization 

a of the Hamiltonian for reasonably large chains (up to approximately 50 spins). 

I These quite precise data allow to carry out a finite size analysis and extrapo- 

i-rt ' late to the thermodynamic limit. On the other hand, 1 — D spin chains have 

5IJ I been analyzed by different techniques, like quantum Monte Carlo methods, cou- 

O ■ pled cluster method (ea;p(S')-method), renormalization group techniques etc., 

i~ I and thus constitute a benchmark problem for comparison of various methods. 

^ ■ From the physics point of view, there are substances occuring in nature having 

KJi I predominantly the structure of a 1 — Z? chain. Examples are CuCh^NC^H^ 

'^ '. or Cu{NH3)iS04,H20 (Mazzi 1955), (Dunitz 1957), (Bernasconi 1981). New 

C^ ' neutron scattering experiments on the quasi 1-D antiferromagnet KCuF^ are 

discussed by Tennant et al. (Tennant 1995). There are recent investigations by 
Dender et al. (Dender 1996) on Copper benzoate also being a linear chain anti- 
ferromagnet. The Heisenberg model has a relation to the Hubbard model. Ac- 
cording to Inui(Inui 1988) the pair hopping term of the Hubbard model can be 
mapped onto the next-to-nearest neighbour interaction of the antiferromagnetic 
Heisenberg (AFH) model. From the theoretical point of view, 1 — D spin chains 
are interesting because its critical properties are predicted by conformal symme- 
try under the hypothesis that this symmetry holds. Finally, 1 — D spin chains in 
the presence of magnetic fields and next-to nearest neighbour coupling display 
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quite a rich phase structure. It is the purpose of this work to consider 1 — D spin 
chains and to present numerical results from exact diagonalization, discuss its 
finite-size analysis and critical behavior. We mainly consider static spin structure 
factors. In particular we analyze chains with next-to-nearest neighbour coupling 
allowing for frustration. 

2 Heisenberg model and its extensions 

Let us start by recalling the definition of various 1 — D spin models considered 
in the following: The classical model, describing next-neighbour interaction is 
given by the Ising model 

N 

H = J^S-S-,^1, (1) 

i=l 

where S^ is the z-component of spin i, taking the values —1, +1. The correspond- 
ing quantum model is the Heisenberg model, 

N 

ff-J^S, -S.+i, (2) 

where S — cr/2 and aa being the Pauli-matrices for the spin 1/2 case. This is 
considered as the standard model of magnetism. In the case J > the model 
describes antiferromagnetism, while J < describes the ferromagnetic system. 
The following generalizations are of physical interest: The XX Z model describes 
the Heisenberg model with an anisotropy in z-direction, 

N 

H^jJ2 StSt^, + SfSf+, + cos{^)StSt+,. (3) 

The Heisenberg model exposed to an external (classical) magnetic field is given 

by 

N N 

// = Ji^S, -S.+i + Jb^B, -S,. (4) 

1=1 i=l 

If one takes into account next-to-nearest neighbour coupling of spins, the model 
becomes 

N N 

iJ = Ji^Si • Si+i + J2^Si • Si+2, a = J2/Ji. (5) 

i=l i=l 

If the sign of Ji and J2 are such that one term favors parallel spins, while 
the other favors antiparallel spins, the system becomes frustrated. A model of 
particular interest with coupling beyond next neighbours is the Haldane-Shastry 
model (Haldane 1988), (Shastry 1988) (with periodic boundary conditions), 

N-l N -. 

H ^ J 2^2^ — ^- — TTTrSm • S„+„. (6) 

■^ "^ sm (niT N) 

n=l 7n=l ^ ' ' 
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Below we will consider combinations of the above cases, e.g., the presence of a 
magnetic field and next-to-nearest neighbour coupling. In the following we will 
consider spin models with periodic boundary conditions. 

2.1 Symmetries and integrability of Heisenberg model 

The isotropic Heisenberg model has a number of symmetries and corresponding 
conserved quantities: There is translational symmetry. The translation operator 
T = exp[ii-'] defines the conserved quantum number of momentum P. There 
is also reflexion symmetry. The total spin S = X]i=i ^i i^ conserved. There 
are more conserved quantities like, e.g., F3, — 2i^„^-^ e-ijkS^^_iSl^S'^j^-^ (Liischer 
1976), (Fabricius 1990). The existence of conserved quantities can be traced 
back to the existence of a set of commuting transfer matrices r(A), satisfying 
[r(A),T(A')] ~ and \T{\),H] = 0. The existence of such transfer matrices for 
the 1 — D Heisenberg model can be proved from inverse scattering theory or via 
the Yang-Baxter equations (Baxter 1980). The set of transfer matrices implies 
the existence of infinitely many conserved quantities like, e.g., the Hamiltonian 

iJ ex 4jlnT(A)L , + const., the total spin squared S^ ex ^wT(A) 



9^ 



A= 



-I a ^ -f));j^^T{X) + const. (Fabricius 1990). Thus the Heisenberg model 

\—i/2 

is an integrable system (Baxter 1980). Analytical solutions have been obtained 
via the Bethe ansatz (Bethe 1931). Using the Bethe ansatz, Hulthen (Hulthen 
1938) has calculated the ground state energy per site in the thermodynamic 
limit eo = 1/4 — In 2 (corresponding to J = 1/2). Analytical information on the 
1 — D model at the critical point can be obtained by making the assumption 
of conformal invariance. Cardy (Cardy 1986), (Cardy 1987) has shown that the 
spectrum of the transfer matrix is determined by conformal invariance. This 
predicts, e.g., the leading term of the finite size behavior of the ground state 
energy, or critical exponents of the spin-spin correlation function. 

Recently, much analytic progress has been made in the Heisenberg and re- 
lated spin models. Romer and Sutherland (Romer 1993) have computed the finite 
size dependence of the correlation function of the AFH model. The dynamical 
correlation function of the Calogero-Sutherland model has been computed by 
Ha (Ha 1994) and the single-particle Green's function of this model has been 
determined by Zirnbauer and Haldane (Zirnbauer 1995). Essler et al. (Essler 
1996) have shown that correlation functions of the XXZ antifcrromagnet near 
the critical magnetic field can be expressed in terms of solutions of Painleve 
differential equations. 



2.2 Quantum numbers of lowest state for finite chains 

The quantum numbers of the lowest lying state for a finite spin chain have been 
analyzed by Fabricius ct al. (Fabricius 1991). For a given total spin S and number 
of spins N one finds the quantum number of momentum P of the lowest lying 
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state to be as follows: For TV even holds 

n/ JO ^ f S*: TV = 4,8,12,- •• ,^, 

For N odd, there is degeneration of momentum 

5(l + l/iV) : A^ = 3,7, II,--- 



P {mod2n)- -'-^sil + im + l-.N^B^Q^U,- 
P (mod 2.) = -. X I ^^^^ _/^/;^^ : TV ^ 5, 9, 13, - - - ^^^ 

The absolute ground state corresponds to 

0: A^ = 4,8,12,--- 



P(^-^)- 1 ^:iV = 2,6,10,■ 



.±^c_l/o^_ /o J (1 + 1/A^) :^^-3,7,ll,- 



^"(^-l/2)-T-/2x|(^l+{/^):^.5,9,13,--- (9) 

These are the quantum numbers corresponding to the isotropic Heisenberg model. 
In the case of next-to-nearest neighbour coupling the situation becomes more 
complicated. Then the quantum numbers of the lowest lying state will depend 
on the model parameters, like a. 



2.3 Magnetization curve and susceptibility 

Magnetization is defined by M =< S > /N, it is a function of the external field 
B. The functional dependence M versus B, the so-called magnetization curve is 
most easily obtained by computing the ground state energy e per spin in a given 
spin sector and using 

de{M) 



dM 



B. (10) 



This follows from 



i/ = i/|B=o+B-S, 

£; = .b|b=o + b-<s>, 

e = ^/7V = e|B=o + B-M. (11) 

The magnetization curve of the AFH with anisotropy at temperature zero is 
shown in Fig.[ 1]. One observes that when magnetization reaches its maximal 
value M = 1/2 (all spins aligned), then the magnetic field B becomes saturated. 
This can be understood as follows: For very long chains, i.e., near the thermody- 
namic limit, and close to saturation (nearly all spins aligned), then the change 
in energy when flipping a single spin becomes independent of the magnetization 
M. Thus B — de/dM — const, and B docs not increase when M increases. The 
magnetization curve for the isotropic AFH model is analytically known to lead- 
ing order in two limiting cases: (a) Weak magnetic field and (b) near saturation. 
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Fig. 1. Magnetization curve computed from Bethe ansatz equations for chain of length 
N=2048. 7 is the anisotropy parameter. Taken from (Karbach 1994a). 



The case of weak magnetic field has been computed by Griffith (Griffith 1964) 
as well as by Yang and Yang (Yang 1966) using the Bethe ansatz and solving 
integral equations numerically. One obtains the expansion 



M{B) ^ - 

IT 



1 ln|lnB| 
1 — fli" a2 9 as 



\tlB h\ B \nB 



(12) 



The lowest order is due to Griffith (Griffith 1964). The coefficient ai — 1/2 has 
been computed by Babujian (Babujian 1983) and a2 = 1/4 has been obtained 
by Lee and Schlottmann (Lee 1987), while a^, is unknown. The magnetization 
curve near saturation is known due to Yang and Yang (Yang 1966), 

V2 



M{B) = --^^B,-B + 0(B, - B), (13) 

and Be ~ 2 for the isotropic case. From the magnetization curve, one obtains 
the magnetic susceptibility x — — an ■ The susceptibility corresponding to the 
magnetization curve of Fig.[ 1] is shown in Fig.[ 2]. Information on antiferro- 
magnetic ordering is encoded in the ground state spin-spin correlation function 
< 0|S(x) • S(0)|0 >. For spin chains being isotropic with respect to coupling of x- 
and y-componcnts of spin, one distinguishes the transverse and the longitudinal 
correlation function, 

ujj{l, TV) == 4 < 0\Sf^-^S{\0 >, j = 1, 2 : transverse, j = 3 : longitudinal. (14) 

Luther and Peschel (Luther 1975) have found the following behavior of the 
correlation functions in the thermodynamic limit 



(-1)' 



COj{l,N)^N^oo^^. (15) 
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Fig. 2. Susceptibility computed from Bethe ansatz equations for chain of length 
N=2048. 7 is the anisotropy parameter. Taken from (Karbach 1994a). 



For the isotropic Heisenberg model the critical exponents r]j obey rjj = 1 for 
j — 1,2,3. This follows from conformal symmetry (Luther 1975), (Bogoliubov 
1986). Also based on conformal symmetry, Affleck et al. (Affleck 1989) have later 
predicted for the isotropic Heisenberg model a logarithmic correction to the long 
range behavior, 

(-1)^ 
< S(x)S{0) >^ A^ ^(Inx)'", cr = l/2. (16) 



3 Finite size scaling 

Critical behavior of a physical theory can exist only in the thermodynamic limit 
(TV ^ oo). Because very few models are soluble in the this limit, it is of great 
practical interest if physical information of the critical system can be obtained 
from a finite number of spins (finite volume). This has been answered by the 
theory of finite size scaling, introduced by Fisher (Fisher 1971). It can be sum- 
marized as follows: Suppose a model becomes critical at a temperature Tc- For 
the reduced temperature t = [T — Tc)/Tc this corresponds to tc = 0. Suppose 
there is an observable, which at the critical point behaves like a power law 

r{t) -^t^o ^, (17) 

where (j) is a. critical exponent. In a finite volume N the observable is r{t,N). 
The theory of finite-size scaling is based on the following assumption relating 
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Fig. 3. Transverse (a) and longitudinal (b) structure factor for N — 4, • • • , 28 and 
different anisotropy parameters 7. Taken from (Karbach 1994a). 



the observable in the thermodynamic limit to the finite volume, 

N 

r{t,N)^N^^rit)fi—), (18) 

where ^(i) is the correlation length in the thermodynamic limit and f{z) is some 
function of a single variable. The correlation length has a critical behavior 

m ^t^o ^. (19) 
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Fig. 4. Longitudinal structure factor versus scahng variable I/3, Eq.( 28), for 
7/7r = 0, 0.1, 0.2, • • • , 0.5 and iV = 4, 6, • • • , 28. 



The finite-size scaling hypothesis can be written in an equivalent form, originally 
proposed by Fisher (Fisher 1971), 

r{t,N) ^N^oo N^^-'giNn ^ ■ (20) 

t=const. 

According to a theorem by Mermin and Wagner (Mermin 1966) the Heisenberg 
model in Z? = 1 and D — 2 dimensions does not have long range order and 
hence does not have a 2nd order phase transition for any finite temperature 
T 7^ 0. The critical exponent v of the correlation length has been computed in 
the thermodynamic limit for the XXZ spin model at the critical point Tc = 
by Suzuki et al. (Suzuki 1992), using the Bethe ansatz. The result is i^ = 1. 



4 Finite size analysis of spin-spin structure factors 

Although the Bethe ansatz in principle allows for an exact calculation of the spin- 
spin correlation function, it is practically very difficult. As an alternative many 
workers have computed numerically correlation functions and structure factors. 
The following standard methods have been applied successfully to spin sys- 
tems: exact diagonalization (Bonner 1964), coupled cluster - or exp(S') - method 
(Bishop 1995), quantum Monte Carlo method (Ceperley 1995), renormalization 
group techniques (White 1993), (Hallberg 1995). In the following we will focus 
on results obtained by exact diagonalization. Due to limitations in computing 
time and storage, exact diagonalization is possible only for small systems (for 
spin 1/2 chain of A^ spins in 1-D, the dimension of Hilbert space is 2^). However, 
in the presence of a strong magnetic field B close to saturation nearly all spins 
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Fig. 5. Finite-size dependence of Aj{^,p, N), Eq.( 31), at momentum p = 7r/2. 



are aligned. Setting up a Hilbert space basis by counting configurations of non- 
aligned (reversed) spins, the dimension will be much smaller. E.g., for A^ = 50 
spins in the sector of magnetization M — 2/5, the dimension of Hilbert space 
is d ~ 42376. Thus spin chains of up to 56 spins have been diagonalized, using 
the Langzos method and binary coding of spin configurations. The advantage 
of exact diagonalization is the high numerical precision allowing to carry out a 
finite size analysis and to extrapolate to the thermodynamic limit. 



Fig. 6. Thermodynamic limit of transverse (T) and longitudinal (L) structure factor: 
Tl, LI finite-size analysis, T2, L2 conjecture of Ref. (Miiller 1981), TS exact result of 
XX-model. (a) 7 = 7r/2, (b) 7 = 0. 
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4.1 Structure factors in momentum space 

Physical information on magnetic ordering is encoded in the spin-spin correla- 
tors. However, from the numerical point of view instead of analyzing the long 
range behavior of the correlation function it is easier to consider its Fourier 
transform, i.e., the spin structure factor in momentum space (Karbach 1993). 
The static spin structure factor is defined by 

Af-l 

S'j(p,Ar) = ^ Wj(;,7V)exp[ip;], p^2Trk/N, fc = 0, • • • , iV - 1. (21) 

(=0 

One has the following correspondence 



x-space: < S{x)S{0) >- 



X 



^tpx\ j_ ^ \ divergent at p = ±7r(22) 

X ^-^ X 

X X 

Thus the long distance behavior of the correlator manifests itself as divergence 
in p-space. Such behavior can be analyzed numerically for finite chains. 



Fig. 7. Scaling function gj in the finite-size scaling ansatz of Eq.( 32) for z — 2. 



4.2 Results for XXZ-model 

For the XXZ-model with anisotropy angle 7, Eq.( 3), the correlation function 
behaves as (Luther 1975) 

a;j(7,0^c(7)^, Vi{l) ^ Va'il) ^ ^ - i/t^, 0<j<tt/2. (23) 
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Note that 771(7) • Vai^f) = 1- The corresponding structure factor behaves as 

S,{j,p)^a,+b,[l-p/7rY>'-\ (24) 

Because 771(7) < 1 and 773(7) > 1, at momentum p ^ n the transverse structure 
factor becomes divergent, while the longitudinal structure factor remains finite. 
This behavior changes when 7^0 (isotropic limit). For finite systems the 
corresponding correlation function and structure factor is given by 

c.,(7,0^5,(7)j^exp[VO(7,A^)], 

S,i^,p) ^ ~a,+b, [{l-p/nf + 1/^1(7,^^)]"^-' , (25) 

where ^j{JtN) is the correlation length in the finite system. The correlation 
length has the following critical behavior when cosh(7') > 1 and 7' ^- (Johnson 
1973) 

1 r 71-2 1 

divergent at 7' = 0. (26) 



C — > o exp 



_27'_ 
The structure factor diverges at p = tt when 7' ^ (Singh 1989) 

S,iy,p^7r,N^^)^y^o^, Ai- 3/2, A3 -2. (27) 

The longitudinal and transverse structure factors for finite systems are shown in 
Fig. [3] 



Fig. 8. Scaling behavior of transverse structure factor at p = tt, compared with fi- 
nite-size scaling ansatz of Eq.( 33). 



Now let us look at the scaling properties of the structure factors. One can ask 
the following questions: (a) Does a scaling variable Lj exist such that Sj (7, p, N) = 
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Sj{'-f, Lj{'y,p, N))7 (b) How docs Sj{j — 0,p,N ~ oo) behave when p ^ it? (c) 
What is the leading finite-size behavior of Sj{'y,p,N) when p < tt? (d) How 
does Sj{'y,p = w, N) behave when N is large and 7^0? Karbach and Miitter 
(Karbach 1993) have suggested to use L^(j,p) as scaling variable, given by 

Lsil^P) = -^TtK fl - (1 -p/TTr'^-"^-'] , P<7r. (28) 

Haldane and Shastry (Haldane 1988), (Shastry 1988) have shown that 

limi3(7,p) = -ln(l-p/7r), p<7r (29) 

7— ►o 

is the exact structure factor of the Haldane-Shastry model, Eq.( 6), with periodic 
boundary conditions. Numerical studies (Karbach 1993), (Karbach 1994b) show 
the following scaling behavior. Fig.[ 4] gives a plot of the longitudinal structure 
factor versus the scaling variable L3. The data points lie close to a line of slope 
one when p < tt. Similar but less impressive scaling holds for the tranverse 
component. Scaling deviations can be further analyzed by looking at 

A(7,P,iV) - Sj{j,p,N)-L,{j,p). (30) 

It is remarkable that at 7 = 7r/2 (XX-model) no finite-size dependence exists. 



Fig. 9. Longitudinal structure factor versus momentum p and magnetization M for 
N = 20,22,---, 28. The ridge occurs at softmode P3(M), Eq.( 34). 



^'3(7 = tt/2,p, N) — 2p/-K. For p < tt one makes the following finite size ansatz 

A,{j,p,N) ^ A,{j,p,N = ^) + ^^^l^, p<n. (31) 
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The finite size behavior is shown in Fig.[ 5] and confirms the leading l/N"^ behav- 
ior when approaching the thermodynamic limit. The case 7 = 7r/2 corresponds 
to the XX-modcl which is exactly soluble. Fig. [6] shows the thermodynamic 
limit of the structure factors. For 7 ~ 7r/2 comparison with the exact result 



Fig. 10. Transverse structure factor at M = 1/4. (a) Finite-size behavior at p = tt, 
Eq.( 35). (b) Momentum dependence for p -^ tt, Eq.( 36). 



shows agreement. The above analysis is not valid when approaching p — tt. The 
following finite-size scaling ansatz can be made in the combined limit N ^ 00, 

P^ TT, 

Sjil^P^N) — 'p^TT^N^oo Sj{-f,p,N = 00) gAl^ z)\z={i-p/^)N=const. (32) 

The scaling function gj is shown in Fig.[ 7]. The transverse structure factor is 
singular ai p = tt. At the singularity the following finite-size ansatz can be made 

5i(7,p = TT, N) -^N^oo n^^ [1 - sin(r7i7r/2)(ai/7V)''i-i] . (33) 

771 - 1 '- 

Numerical data compared with this ansatz are shown in Fig.[ 8]. 

5 AFH model with next-to-nearest neighbour couphng 
and external magnetic field 

Antiferromagnetic order is manifested in the Heisenberg model via singularities 
of structure factors when p — > tt, iV ^ 00 and z = {1 — p/tt)N = const. Then 
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the transverse structure factor becomes infinite, while the longitudinal structure 
factor remains finite. This antiferromagnetic order can be destructed (a) by 
switching on an external uniform magnetic field B or (b) by frustration (i.e., 
nearest neighbour and next-to-nearest neighbour coupling do not favour the 
same alignment pattern). Such will be examined in this section. As a result it 
turns out that there are similarities between the anisotropic chain (XXZ) and 
the chain with next-to-nearest neighbour coupling. 



Fig. 11. Longitudinal structure factor at M = 1/4. (a) Finite-size behavior at softmode 
P3(M), Eq.( 38). (b) Momentum dependence for p -^ pa.{M), Eq.( 39). 



5.1 Presence of external magnetic field 

When considering the Heisenberg Hamiltonian in the presence of a uniform ex- 
ternal magnetic field B, Eq.( 4), (with B^ = B = Be^ and Jb = 1), one observes 
in the spectrum of the Hamiltonian the following property: Zero-frequency modes 
^{Psoft) = emerge at non-zero "soft-mode" momenta Pso/t • These modes show 
up in the dynamical stucture factors Si{uj,p) and S'3(cj,p) (Fledderjohann 1996). 
The momenta depend on the field B, or the magnetization M, respectively. One 
has 

Pi{M) = 27rM, paiM) = n{l - 2M). (34) 

The critical exponents rjj are supposed to govern the long distance behavior of 
the static structure factors and the infra-red behavior of the dynamical struc- 
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ture factors. Now rji — rji{M) become field dependent. Provided that the low- 
lying excitation spectrum is governed by conformal symmetry, one can make 
analytical predictions. Numerical studies (Karbach 1995), (Fledderjohann 1996) 
show the following behavior. In Fig.[ 9] the longitudinal structure factor is plot- 
ted versus momentum and magnetization. The results correspond to chains of 
A'' = 20, 22, • • • , 28 spins. The structure factor is a smooth function except for a 
ridge at the soft-mode. One observes that the soft-mode moves with magnetiza- 
tion. On the other hand, the transverse structure factor diverges at p — n when 
N ^ oo like 

Si{p^7:,M,N)^N^^N'-^'^^\ (35) 

For M = 1/4 one finds, e.g., rji — 0.65, shown in Fig.[ 10a]. The exponent rji 
governs also the approach to the singularity at p = tt, 

Slip, M, TV - oo) -^p^^ (1 - p/tt)"!'^)-!, (36) 

which is shown in Fig[ 10b] . When approaching the soft-mode momentum pi (M), 
one finds the behavior 

5i(p,M,iV-^)^p^p,(M)±o (l-p/pi(M))''*('^^-\ (37) 

Numerical results indicate that ?7f (M) and ri^(M) do not necessarily agree. 
The corresponding behavior of the longitudinal stucture factor is as follows. The 
finite size behavior at the soft-mode momentum p^lM) is given by 

S3ip3iM),M,N)^N^ooN^-^^'^^^. (38) 

This is displayed in Fig.[ 11a]. Approaching the soft-mode from the left one has 



Fig. 12. Magnetization curve for a = and a — 1/4. The solid line is the Bethe ansatz 
solution with N=2048 spins. 
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Ssip, M,N = ^) ^p^p3(Af)-o (1 - p/p3(M))''-^('^^-^ , (39) 

shown in Fig.[ lib]. There is indication for rj^{M) — rj^{M). However, if r7;^(M) 
and f]^ (M) agree is an open question. 

5.2 Next-to-nearest neighbour coupling 

Now let us consider the model Hamiltonian including next-to-nearest neigh- 
bour coupling and the presence of a magnetic field, given by a combination of 
Eqs.( 5, 4). This model has a rich phase structure, depending on the parame- 
ters Ji, J2 and B. There is a ferromagnetic phase, an antiferromagnetic phase 
and a frustrated phase (Farnell 1994). Haldane (Haldane 1982) has shown that 
within the antiferromagnetic phase there is a phase transition at ac between a 
spin-fluid phase (a < ac) and a dinier phase {a > ac). The transition point has 
been determined to be ac = 0.24 (Okamoto 1992). Related to this is Gluzman's 
(Gluzman 1994) prediction of the occurence of a plateau in the magnetization 
curve. Moreover, it is remarkable that the model is exactly soluble at a = 1/2, 
B = 0, called Majumdar-Ghosh point (Majumdar 1969). One can ask: Where 
are the observable differences in the model with and without next-to-nearest 
neighbour coupling, i.e., a 7^ and a = 07 Numerical studies by Schmidt et 
al. (Schmidt 1996) and Gerhardt et al. (Gerhardt 1996) give the following in- 
formation. Firstly, one observes in the magnetization curve near saturation the 
following characteristic behavior, 

M{a = 1/4, B) — >B^B., 1/2 - const, x (B, - B)^/^, (40) 

in contrast to the square root behavior for a = 0, given by Eq.( 13). This is shown 
in Fig.[ 12]. In the presence of an external magnetic field and next-to-nearest 
neighbour coupling, the structure factors develop singularities at soft-mode mo- 
menta, given by Eq.( 34). Now the strength of the singularity will depend on the 
the frustration parameter a. The transverse structure factor and the longitudi- 
nal structure factor for fixed values of a but different values of magnetization 
M are shown in Fig.[ 13a] and [ 13b], respectively. Fig.[ 14] displays the struc- 
ture factors for fixed M but different values of a. The critical exponents rji will 
now depend on M and a. The finite-size dependence of the transverse structure 
factor at p = TT is described by 

Si{a,p ^ TT, M, N) — >N^oo Ai{a, M) + Bi{a, M)N^-'^'^"'^\ (41) 

E.g., one finds T]i{a ^0,M ^ 1/4) == 0.65, while f]i{a = 1/4, M == 1/4) ^ 1.16. 
Similarly, the longitudinal structure factor at the soft-mode p3(M) behaves as 

S3{a,p = P3{M),M, N) -^N^oo A^ia, M) + B^{a, M)N^-'^^'^''-^\ (42) 

E.g., one finds 7/3(0; = 0,Af = 1/4) = 1.5, while 773(0: == 1/4, M == 1/4) == 0.84. 
Those critical exponents extracted from Eqs.( 41, 42) are given in Tab.[ 1] as 
a function of a. Based on the assumption that conformal symmetry holds, one 
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Fig. 13. Transverse structure factor (a) and longitudinal structure factor (b) versus 
momentum p for fixed magnetization M = 1/4 but different values of a. 



would expect '7i • ??3 = 1 • The numerical data for most values of a deviate from 
the expected result by a few percent. However, special attention should be paid 
to the region a ^ —0.2 and the region a > 0.4. For a = —0.2 Eq.( 42) does 
not yield a value of 7^3. It turns out that the finite-size depence is very weak 
in this region (B^i^a, M) is very small) and may even vanish for some value of 
a. Secondly, for a > 0.4, rji ■ rj^ deviates substantially from the value one. The 
reason is most likely the following: The structure factors have been computed in 
the sector corresponding to quantum numbers of the Heisenberg model, Eq.( 9). 
For a > 0.4, those are presumably no longer the quantum numbers of the ground 
state, hence the critical exponents rji do not correspond to the ground state but 
to some excited state. 
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Fig. 14. Transverse structure factor (a) and longitudinal structure factor (b) versus 
momentum p for fixed a = 1/4 but different values of M. 



Table 1. Critical exponents rji for M = 1/4 versus a 



a 


-0.5 


-0.4 


-0.3 


-0.2 


-0.1 





0.1 


0.2 


0.25 


0.3 


0.4 


0.5 


m 


0.39 


0.42 


0.45 


0.50 


0.56 


0.65 


0.78 


1.00 


1.16 


1.36 


2.74 


3.53 


V3 


2.52 


2.45 


2.34 


? 


1.66 


1.49 


1.26 


0.98 


0.84 


0.70 


0.47 


0.32 


Vi -Vi 


0.98 


1.03 


1.06 


? 


0.93 


0.97 


0.98 


0.98 


0.97 


0.94 


1.28 


1.13 
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6 Summary and outlook 

The 1-D spin 1/2 AFH model in the presence of an external magnetic field and 
with next-to- nearest neighbour coupling can be solved " exactly" by diagonaliza- 
tion (via Langzos) of chains of up to 50 spins (depending on M). For the ground 
state properties of the model this yields data with high numerical precision al- 
lowing for finite-size scaling analysis and extraction of critical exponents. Here 
we have concentrated on static spin structure factors. One finds a different be- 
havior for a < ac and a > ac- The difference is observed in the strength of the 
singularities of structure factors and in the magnetization curve near saturation. 
Much interesting physics is to be explored when going into D=2,3 dimensions 
(e.g., prediction of plateaus in the magnetization curve). Exact diagonalization 
is presently limited in D=2 to 7 x 7 spins (with helical boundary conditions). 
To analyze critical behavior in D=2, new methods yielding numerically precise 
results for a larger number of spins are needed. 
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